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Abstract 

The interaction of metal-molecule-metal junctions with light is considered within a simple generic model. 
We show, for the first time, that light induced current in unbiased junctions can take place when the bridging 
molecule is characterized by a strong charge-transfer transition. The same model shows current induced 
light emission under potential bias that exceeds the molecular excitation energy. Results based on realistic 
estimates of molecular-lead coupling and molecule-radiation field interaction suggest that both effects should 
be observable. 

PACS numbers: 73.23.-b,78.20.Bh,78.20. Jq,78.60.Fi,78.67.-n 
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Molecular conduction nan o-j unctions have been under intense study for some time. A class 
of molecules not yet investigated in this context are those characterized by strong charge-transfer 
transitions into their first excited state. The dipole moment of such molecules changes consid- 
erably in such transitions [lj, expressing a strong shift of the electronic charge distribution. In 
the independent electron picture, the charge transfer nature of the first excited state implies that 
either the highest occupied, or the lowest unoccupied, molecular orbitals (LUMO and HOMO) is 
dominated by atomic orbitals of larger amplitude (and better overlap with metal orbitals) on one 
side of the molecule than on the other and therefore stronger coupling to one of the leads. We 
show that when such molecular wire connects between two metal leads with the charge transfer 
direction approximately parallel to its axis, optical pumping may create an internal driving force 
for charge flow between the leads. Such optical-resonance activation of current flow is different from 
previously demonstrated adiabatic pumpingQ and from previously discussed strong field optical 
control mechanisms j3, 0]. 

The opposite phenomenon: light emission in current carrying tunnel junctions has already been 
demonstrated. Emission from bare junction : ;5|| was attributed to surface plasmonsyl. Emission 

PQ 

from molecular nanojunctions[7], la] was n °t so f ar discussed theoretically. 

Our model provides a simple framework for treating both light induced current and current 
induced light in molecular junctions. Following its introduction below we present results of self 
consistent numerical solutions as well as approximate expressions for the corresponding electron and 
photon fluxes. We examine, using reasonable experimentally-based parameters, the magnitudes of 
these effects and conclude that both phenomena can be observable under fairly general conditions. 

In our model (Fig. ^) the molecule is represented by its HOMO, |1 >, and LUMO, |2 >, with 
energies E\ and £2 and gap £21 = £2 — £1, positioned between two leads represented by free electron 
reservoirs L and R characterized by the electronic chemical potentials hl and fiR. fii — fJ-R = 6*3? is 
the voltage bias. In the independent electron picture molecular excitation corresponds to transfer 
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FIG. 1: The molecular junction model. The right (R — {\r >}) and left [L = {\l >}) manifolds represent 
the leads' electronic states; fiL and fin are the corresponding electrochemical potentials. The molecule is 
represented by its HOMO and LUMO levels. 

of an electron between levels |1 > and |2 >. The corresponding Hamiltonian is 

H = H + V M + V n + Vp = Hj + V n + V p (1) 
#o = E e mclnCm+ E e k c\c k + h E uj a ala a (2) 

m=l,2 k£{L,R} a 

Vm= E E (^4^ + H.c.) (3) 

K=L,Rm=l,2;keK 

Vn= E E ( W*)4^k + Hx -) (4) 

K=L,Rk^k>eK 



V^^aoCgCi + H.c.) case a 

So ( V"i P) a Q 4ci + H.c. ) case b 



(5) 



where H.c. denotes Hermitian conjugate. Hq in © contains additively terms that correspond 
to the isolated molecule, the free leads and the radiation field. Eqs. ©-© describe the coupling 
between these subsystems. The operators c (c') and a (a)) are annihilation (creation) operators for 
electrons and photons. Vm is the standard electron transfer coupling that gives rise to net current 
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in the biased junction and Vjy describes energy transfer between the molecule and electron-hole 
excitations in the leads. The latter interaction strongly affects the lifetime of excited molecules near 
metal surfaces and should be inlcuded in any analysis of phenomenainvolving resonance radiation 
and molecular conduction junctions. The molecule-radiation field coupling Vp, Eq. ©, will be 
taken in two forms: The form describes driving of the junction by the electromagnetic field 
mode a = 0, Eq. (JSJa) is used when addressing the issue of spontaneous light emission from 
current carrying junctions. We limit ourselves to near resonance processes pertaining to linear 
spectroscopy. This justifies the use of the rotating wave approximation (RWA) in Eq. (JSJ). Note 
that the radiative coupling coefficients reflect properties of the local electromagnetic field at the 
molecular bridge which depend on the metallic boundary conditions. 

In the Keldysh non-equilibrium Green function approach the steady state currents through the 
bridge are obtained from 



were the Green functions (GFs) G and the self energies (SEs) £ are defined in the bridge subspace. 
The subscript B corresponds to a particular relaxation processes represented by the SE S^^. At 
steady state the absorption and emission photon fluxes, I a & s and I em , the non radiative relaxation 
flux, In, and the source-drain current, I st j, come into balance. The SEs needed for their evalua- 
tion are calculated as sums of independent contributions associated with the different relaxation 
processes (the non-crossing approximation) 



where Y, x (X = ML, MR, P, NL, NR) is the SE associated with the coupling VW in Eqs. 
On the Keldysh contour these SEs are 2x2 matrices in the bridge space 





(7) 




(8) 





(9) 



a 



F a (T 2 ,Tl)G 2 2(Tl,T 2 ) 







(10) 



X 
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here and below K = L,R denotes the left and right heads, m,m' = 1,2 and in = 25 m ^ + 
&m,2- 9k an d F a are free electron and photon GFs in state k and mode a, respectively. The 
retarded/advanced, lesser and greater projections of these SEs on the real time axis are obtained 
using the Langreth formulas and can be expressed at steady state situations in energy space. In 
the wide-band approximation and assuming that the HOMO and LUMO are not mixed by their 
interactions with the leads, the SEs associated with electron exchange between molecule and leads 
have the familiar forms 



^ r MK,mm> ~ [^MK,m'm] ~ — ^mm'T W,m/2 
S MJ(,ram' = idmm'fK(E)T M K,m 

—id mm i 

[1 - fK{E)]Tj\4K 



keK 

f K (E) = [exp {(E - fi K )/k B T} + l]" 1 



(11a) 
(lib) 
(11c) 
(lid) 

(lie) 



The SEs associated with the non-radiative and radiative energy relaxation processes are similarly 
obtained in the forms 



/dto 
■^-B nk (uj, fi K )G^ m (E + w) 

/d 
—B NK (u, fi^G^iE - w) 

— C NK (E,u)f K (E)[l - f K (E + co)] 



(NK) 



C NK (E,u) = (2ir) 2 E \ V k« 

k^k'tK 

x 5{E - e k )S(E + u - e k .) 



(12a) 
(12b) 
(12c) 



(12d) 



and 



a 

(l + N a )G<(E + u 



x 



a 

N a G> 2 {E + uo a ) 



x 





N a G<(E-u c 





(1 + NJG^E - u a ) 



(13a) 



(13b) 
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where N a is the number of photons in mode a. The corresponding retarded and advanced SEs are 
more difficult to calculate from the Langreth rules. For simplicity we assume, in the spirit of the 
wide band approximation, that all diagonal components of £ are purely imaginary Consequently 

[£>(£)-£<(£)] 



£ r (£) = [£ a (£)] T 



-Ir 
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(14) 



We use (fTl)) to get the retarded and advanced GFs according to G 7 mm ,(E) = (G^, m (£'))* = 
[E — e m — SJ^ m (i?)] _1 S mm r. G <,y needed in © are then obtained from the Keldysh equation 
G <>> (E) = G r (E)Yl < ' > (E)G a (E). After evaluating these matrices we can use Eq. © with the 
appropriate self energy to calculate the desired flux. We are particularly interested in I em and I S( i- 
With regard to the radiative interactions we consider two models: 

(A) To study the effect of pumping the junction with light we use a model with a single photon 

mode of frequency loq and disregard all other modes, as described by Eq. (j^H)- In the semiclassical 

* (P) 

limit for the radiation field and in the RWA we can set Nq = 1 and take the coupling Vg as 
a product of the local electric field amplitude and the molecular transition dipole 9] . The charge 
transfer transition on the bridge is expressed in this model by taking LUMO bridge level more 
strongly coupled to one lead than to the other. For the HOMO this inequality is assumed weaker 
or reversed. The flux of interest is I s d under illumination. 

(B) To describe current driven spontaneous emission we use the interaction (0)) with all radiation 
field modes taken in their vacuum state. The frequency resolved emission, I' em {uj) = dl em {ui) / duj 
is obtained from (jHJ) using the frequency resolved self energies obtained from (|13|) 



7pM 
2irpp(io) 

7pM 



G<(E + u) 







Gf^E-to) 



2irpp(uj) 
with the radiative width 

7P (w) = 2vr ]T Wp 5(oj - oj a ) = 2vr VP ' Pp(lu) 



(15a) 
(15b) 



(16) 



and the density of photon modes, Pp(uj) = lo 2 /(tt 2 c 3 ) where c is the speed of light. The total 
emission flux is found from i^n = Jo°° ^ C^)' 

In general, the coupled equations for the self energies and the Green functions have to be solved 
self consistently. The numerical results displayed below were obtained from such an iterative solu- 
tion. More transparent closed forms may be written under the additional simplifying assumptions: 



(a) £21 is large relative to the total widths of levels 1 and 2; (b) Cnk(E,u>); K = L, R, Eq. (|12d|) . 
is taken constant independent of E and u; (c) 7p(w), Eq. (|16|). is also assumed constant and (d) 
keep only teh lowest (second) order in the radiative coupling Vp. It is then possible to express 
the relevant self energies and the resulting currents in terms of the electronic populations n\ and 
n 2 of levels 1 and 2 using n m = (2iri)~ 1 /_ °° dEG^ im (E), m = 1,2. In addition, for case (A) we 
are interested in the small or no voltage regime, whereupon to a good approximation m = 1 and 
n 2 = 0. This leads, for model (A), to 



Isd = I 5^ r ML,mG'^ m (£ , )rMij,mG^ m (^) [Jl(E) - fp(E)] (17) 

/ 



m=l,2 

(P) |2 /* T MLil T MR>2 f L (E - w )[l - - T ML)2 T MRil f R (E - uj )[1 - f L (E)} 



2nh [(E - e 2 y + (r 2 /2)2] x [(S - wq - ei) 2 + (rx/2) 2 ] 

The first term on the right is the usual Landauer term that vanishes at zero potential bias, fi = fp. 
The second term shows explicitly the effect of pumping. In the absence of bias we get 



sd 



y(P) 2 



dE f(E - u )[l - f(E)] {Tml,iT M r,2 - Tml^mra} , , 

2^n[{E-e 2 f + {T 2 /2f]y.[{E-u -e l f + (T l /2f] [ > 

where / = /l = fp- For our choice of parameters below the analytical results 1)17(1 and (|18|) provide 
excellent approximations to the full self-consistent calculation. Furthermore, when u is not too far 
from its resonance value £21 we can approximate f(E — coq)[1 — f(E)] by 1 to get 

2 

x ! ~F~F ' ( 19 ) 



isd 



y(P) 



n (e 2 -w -ei) 2 + (r/2) 2 r!r 2 

Eqs. (|18[) and (|19|l show explicitly how asymmetry in the HOMO and LUMO couplings to the 
metal electrodes leads to photocurrent in the present model. 

For model (B) the frequency resolved emission is obtained from © and 1)15)1 in the form 



I'em(u) 



7P (w) f + °° dE f f L (E + u;)T MLt2 + fp(E + Lo)T MRt2 [1 - f L (E)]T M L,i + [1 - f R (E)]T MR 



x 



h 2n { {E + u-e 2 y + {T 2 /2y (E - £l ) 2 + (Fx/2) 2 

(20) 



and the integrated emission is obtained, using the same approximations as before, in the anticipated 
form 

jtot _ 7p(£2i) n s / 91 x 

iem, - * n 2 (l-nij [21) 

The level populations n\ and n 2 should be obtained from the full self consistent calculation. For 
$ > £21 both LUMO and HOMO bridge levels are well inside the energy window between the 
leads' chemical potentials. A good approximation for these populations is (written for the case of 



negatively biased left electrode) ni = Tml, 2/T2; n\ = Tml,i/^i + ^{Bn + 7p)/Ti- (BP then 
leads to 

rtot _ IP ^ML,2^MR,1 fc)t) s 

lem ~ h 1^2 [ZZ) 

Comparing to (|18|) we see that tailoring the molecule-leads coupling asymmetry such that resonance 
radiation will induce electron flow in a given direction enhances light emission under such bias that 
leads to electron current in the opposite direction. 
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FIG. 2: The light induced current in the model of Fig.Q]in the absence of external potential bias. Full (red) 
line: T = 300 K, e 2 i = 2 eV, Y M la = T M r,i = 0.2 eV, T M l,2 = 0.02 eV, T M r,2 = 0.3 eV, 7P = lO" 6 eV, 
Bnl — — 0.1 eV and V^ P ' ) = 10 -3 eV. Dashed (blue) line — same parameters except TMK.m X 3 

(m = 1, 2) is used. 

Numerical results obtained from the self consistent calculation are shown in Figs. |2 and |31 

Figure 121 shows the light induced current obtained using a set of reasonable junction parameters. 

(P) — ,: i 

In particular, for a molecule with transition dipole moment ~ ID the choice Vg = 10 eV 
corresponds to an incident field intensity ~ 10 8 watt/cm 2 . 
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FIG. 3: The integrated photon emission rate (full line, red) and yield (dashed line, blue) plotted against the 
bias voltage. Parameters are as in Fig. ^except TuK,m = 0.1 eV, K — L, R, m = 1, 2. 

Fig. |21 shows the integrated emission as well as the predicted photon emission yield. The later 
is defined as the ratio between the emitted photon current and the electronic current passed by 
the junction — both displayed as functions of the bias voltage. 

These results, obtained under rather conservative choices of damping parameters, indicate that 
both light induced current and the molecular mechanism (as opposed to the surface plasmon 
mechanism) for current induced emission in tunnel junctions can lead to measurable signals. This 
molecular mechanism can in principle be distinguished from the plasmon mechanism for light 
emission by the frequency dependence of the resolved emission, as implied by Eq. q2UJI . 

In conclusion, we have described a model that accounts for observed current induced light 
emission from molecular tunnel junctions and provides the tools for determining the intensity and 
yield of such emission as functions of key junction parameters. The same model predict that 
resonant light induced current can occur in junctions employing molecular bridges with strong 
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charge-transfer transition to their first excited state. We have concluded, using reasonable param- 
eter, that light driven electronic currents and current driven light emission are realistic possibilities. 
The linewidths and lineshapes associated with these signal can also be analyzed by the same the- 
oretical framework. Correlating observations with predictions made in this paper should help the 
interpretation of future experiments in this direction. 
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